n=2000; %总的分子数
B=ones(1,n);    %起初所有分子在左边，其值为1
for j=1:10000   %经过10000次的随机选择
    k=ceil(rand(1)*n);  %每次是第n个分子穿过隔板
    B(k)=-B(k); %穿过隔板后改变符号
    cc=length(find(B==1));  %统计左边的分子数
    c(j)=cc;
end
subplot(2,1,1)
plot(c)
subplot(2,1,2)
x=0:10000;
y=n*0.5*(1+exp(-2*x./n));   %用解析式计算右边分子数
plot(x)